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The possibility of appearance of spin polarized states in strongly asymmetric nuclear matter is 
analyzed within the framework of a Fermi liquid theory with the Skyrme effective interaction. The 
zero temperature dependence of the neutron and proton spin polarization parameters as functions 
of density is found for SLy4 and SLy5 effective forces, ft is shown that at some critical density 
strongly asymmetric nuclear matter undergoes a phase transition to the state with the oppositely 
directed spins of neutrons and protons while the state with the same direction of spins does not 
appear. In comparison with neutron matter, even small admixture of protons strongly decreases 
the threshold density of spin instability. It is clarified that protons become totally polarized within 
very narrow density domain while the density profile of the neutron spin polarization parameter is 
characterized by the appearance of long tails near the transition density. 

PACS numbers: 21.65.+f; 75.25.+z; 71.10.Ay 



I. INTRODUCTION 

The spontaneous appearance of spin polarized states 
in nuclear matter is the topic of a great current interest 
due to relevance in astrophysics. In particular, the effects 
of spin correlations in the medium strongly influence the 
neutrino cross section and neutrino luminosity. Hence, 
depending on whether nuclear matter is spin polarized 
or not, drastically different scenarios of supernova explo- 
sion and cooling of neutron stars can be realized. An- 
other aspect relates to pulsars, which are considered to 
be rapidly rotating neutron stars, surrounded by strong 
magnetic field. There is still no general consensus regard- 
ing the mechanism to generate such strong magnetic field 
of a neutron star. One of the hypotheses is that magnetic 
field can be produced by a spontaneous ordering of spins 
in the dense stellar core. 

The possibility of a phase transition of normal nuclear 
matter to the ferromagnetic state was studied by many 
authors. In the gas model of hard spheres, neutron mat- 
ter becomes ferromagnetic at g « 0.41 fm~^ (J. It was 
found in Refs. 0, Q that the inclusion of long-range at- 
traction significantly increases the ferromagnetic transi- 
tion density (e.g., up to g « 2.3 fm~'^ in the Brueck- 
ner theory with a simple central potential and hard core 
only for singlet spin states Q)- By determining mag- 
netic susceptibility with Skyrme effective forces, it was 
shown in Ref. |^ that the ferromagnetic transition oc- 
curs at g « 0.18-0.26 fm~'^. The Fermi liquid crite- 
rion for the ferromagnetic instability in neutron matter 
with the Skyrme interaction is reached ai g w 2~4f3o 0, 
where qq = 0.16 fm^^ is the nuclear matter saturation 



density. The general conditions on the parameters of 
neutron-neutron interaction, which result in a magnet- 
ically ordered state of neutron matter, were formulated 
in Ref. 6]. Spin correlations in dense neutron matter 
were studied within the relativistic Dirac-Hartree-Fock 
approach with the effective nucleon-meson Lagrangian in 
Ref. 0, predicting the ferromagnetic transition at sev- 
eral times nuclear matter saturation density. The impor- 
tance of the Fock exchange term in the relativistic mean- 
field approach for the occurrence of ferromagnetism in 
nuclear matter was established in Ref. 1^. The stabil- 
ity of strongly asymmetric nuclear matter with respect 
to spin fluctuations was investigated in Ref. ^ , where it 
was shown that the system with localized protons can de- 
velop a spontaneous polarization, if the neutron-proton 
spin interaction exceeds some threshold value. This con- 
clusion was confirmed also by calculations within the rel- 
ativistic Dirac-Hartree-Fock approach to strongly asym- 
metric nuclear matter |1C| . Competition between ferro- 
magnetic (FM) and antiferromagnetic (AFM) spin or- 
dering in symmetric nuclear matter with the Skyrme ef- 
fective interaction was studied in Ref. [Tlj . where it was 
clarified that FM spin state is thermodynamically prefer- 
able over AFM one for all relevant densities. 

For the models with realistic nucleon-nucleon (NN) in- 
teraction, the ferromagnetic phase transition seems to be 
suppressed up to densities well above qq In 
particular, no evidence of ferromagnetic instability has 
been found in recent studies of neutron matter [l5j and 
asymmetric nuclear matter 16] within the Brueckner- 
Hartree-Fock approximation with realistic Nijmegen II, 
Reid93 and Nijmegen NSC97e NN interactions. The 
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same conclusion was obtained in Ref. [17|, where mag- 
netic susceptibility of neutron matter was calculated with 
the use of the Argonne vis two-body potential and Ur- 
bana IX three-body potential. 

Here we continue the study of spin polarizability of 
nuclear matter with the use of an effective NN interac- 
tion. As a framework of consideration, a Fermi liquid 
(FL) description of nuclear matter is chosen As 
a potential of NN interaction, we use the Skyrme effec- 
tive interaction, utilized earlier in a number of contexts 
for nuclear matter calculations [23-|2^. The main em- 
phasis will be laid on strongly asymmetric nuclear mat- 
ter and neutron matter as its limiting case. We explore 
the possibility of FM and AFM phase transitions in nu- 
clear matter, when the spins of protons and neutrons are 
aligned in the same direction or in the opposite direc- 
tion, respectively. In contrast to the approach, based on 
the calculation of magnetic susceptibility, we obtain the 
self-consistent equations for the FM and AFM spin order 
parameters and find their solutions at zero temperature. 
This allows us to determine not only the critical density 
of instability with respect to spin fluctuations, but to es- 
tablish the density dependence of the order parameters 
and to clarify the question of thermodynamic stability of 
various phases. 

Note that we consider the thermodynamic properties 
of spin polarized states in nuclear matter up to the high 
density region relevant for astrophysics. Nevertheless, we 
take into account the nucleon degrees of freedom only, 
although other degrees of freedom, such as pions, hyper- 
ons, kaons, or quarks could be important at such high 
densities. 



II. BASIC EQUATIONS 

The normal states of nuclear matter are described by 
the normal distribution function of nucleons /kik2 — 
Tr ga'^^a^i, where k = (p, fi, r), p is momentum, (t(t) 
is the projection of spin (isospin) on the third axis, and 
g is the density matrix of the system. The energy of 
the system is specified as a functional of the distribution 
function f,E = E{f ), and determines the single particle 
energy 

eKiK2(/) = 7T7 • (1) 

The self-consistent matrix equation for determining the 
distribution function / follows from the minimum condi- 
tion of the thermodynamic potential 18] and is 

/ = {exp(ro£ + + 1}"' = {exp(roO + 1}"' • (2) 

Here the quantities s and I4 are matrices in the space 
of K variables, with l4ftiK2 = ^4rii5fti«;2 ("^i = ":-P): 
Yo = 1/T, Kin = and Y^p = ^fi'^p/T the La- 

grange multipliers, /ij^ and /ip the chemical potentials 
of neutrons and protons, and T the temperature. Since 



it is assumed to consider a nuclear system with an ex- 
cess of neutrons, the positive isospin projection is as- 
signed to neutrons. This is different from the formalism 
of Ref. [TH , aimed to investigate symmetric nuclear mat- 
ter. Further we shall study the possibility of formation 
of various types of spin ordering (ferromagnetic and an- 
tiferromagnetic) in nuclear matter. 

The normal distribution function can be expanded in 
the Pauli matrices at and in spin and isospin spaces 

/(P) = /oo{p)CTo'ro + /3o(p)CT3'ro (3) 

+ /03(P)(^0T3 + /33(P)(^3T3- 

For the energy functional invariant with respect to ro- 
tations in spin and isospin spaces, the structure of the 
single particle energy is similar to the structure of the 
distribution function /: 

e(p) = eoo(p)CToro + e3o(p)CT3To (4) 

+ e03(p)'70T3 + e33{p)(J3T3. 

Using Eqs. Q, O, one can express evidently the distri- 
bution functions /qq , /30 , /03 , /33 in terms of the quanti- 
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Here n{uj) ~ {exp(yo'.^) + 1}""'^ and 

^+,+ = Coo + £,30 + £03 + £33, 
^+,- = Coo + £,30 — £03 — £33, 
= £00 — £30 + £03 — £33, 
<^-,- = £00 — £30 — £03 + £33, 

where 

£00 = £00 — A^oo) £30 = £30, 
£03 — £03 — Mo3' £33 — £33, 

_ r n ^ h^p _ A*ri H-p 

A^OO — 2 ' ~ 2 ' 

As follows from the structure of the distribution functions 
/, the quantity uj±^±, being the exponent in the Fermi 
distribution function n, plays the role of the quasiparticle 
spectrum. We consider the case when the spectrum is 
four-fold split due to the spin and isospin dependence of 
the single particle energy £(p) in Eq. Q. The branches 
^±,+ correspond to neutrons with spin up and spin down, 
and the branches to±.- to protons with spin up and spin 
down. 
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The distribution functions / should 
Uzation conditions 


satisfy the norma- 


£00 (p) = 
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^ X! /o3(p) = Q71 - Qp = ag, 
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£03 (p) = 


^^C/2(k)/o3(q), 
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-yZ^faoip) = - Qi= Apii, 
p 


(8) 


£33 (p) = 


^5]f/3(k)/33(q). 



:^2^/33(p) = (ftiT + (?pi) - (ftU + = ^QU- (9) 
p 

Here a is the isospin asymmetry parameter, gi„|, p„| and 
PpT I ^?pi ^'^'^ ^^'^ neutron and proton number densities 
with spin up and spin down, respectively; gi = gn'\ + gp| 
and gi = gni + gpi are the nucleon densities with spin 
up and spin down. The quantities Ap|-f and may 
be regarded as FM and AFM spin order parameters. In- 
deed, in symmetric nuclear matter, if all nucleon spins 
are aligned in one direction (totally polarized FM spin 
state), then Ag|| — g and Agi|| — 0; if all neutron spins 
are aligned in one direction and all proton spins in the 
opposite one (totally polarized AFM spin state), then 
Agi|x — g and Apu = 0. In turn, from Eqs. one 
can find the neutron and proton number densities with 
spin up and spin down as functions of the total density 
p, isospin excess 5g = ag, and FM and AFM order pa- 
rameters A(3|| and Ag||: 
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In order to characterize spin ordering in the neutron 
and proton subsystems, it is convenient to introduce neu- 
tron and proton spin polarization parameters 



g-n 



n„ 



(10) 



The expressions for the spin order parameters Apjf and 
A^j|| through the spin polarization parameters read 

= gnn„ + gpUp, Ag^i = gi„n„ - gpiip. 

To obtain the self-consistent equations, we specify the 
energy functional of the system in the form 

Eif) ^ Eoif) + E,r.t{f), (11) 
Eoif) = 4^eo(p)/oo(p), £o(p) = 



2mn 



E.Mif) = 2^{eoo(p)/oo(p) +£3o(p)/3o(p) 



Here niQ is the bare mass of a nucleon, ?7o(k), C^3(k) 
are the normal FL amplitudes, and £ooj £30j £o3j £33 are 
the FL corrections to the free single particle spectrum. 
Further we do not take into account the effective tensor 
forces, which lead to coupling of the momentum and spin 
degrees of freedom [Ij, |25|, 12^ , and, correspondingly, to 
anisotropy in the momentum dependence of FL ampli- 
tudes with respect to the spin polarization axis. Using 
Eqs. 10 and (|ll|l . we get the self-consistent equations in 
the form 

Coo(p) = £o(p) + £oo(p) - Moo- C3o(p) = £3o(p), (12) 
C03(P) = £03 (P) - ^J■03, 63 (p) = £33 (p)- 

To obtain numerical results, we use the Skyrme effective 
interaction. In the case of Skyrme forces the normal FL 
amplitudes read 

C/o(k) = 6to + tsg" + ^[3ti + t2(5 + 4a;2)]k2, (13) 

f/i(k) = -2io(l - 2xo) ~ \t3g\\ - 2x3) 



- - 2xi) - t2(l + 2x2)]k^ = a 
C/2(k) = -2io(l + 2xo) - + 2x3) 

- -|[ii(l + 2xi) -t2(l + 2x2)]k2. 



6k^ 



t/3(k) = -2to - ^<3g^ - ^(ii -<2)k2 = c + dk^ 



where ti,Xi,l3 are the phenomenological constants, char- 
acterizing a given parameterization of the Skyrme forces. 
In the numerical calculations we shall use SLy4 and SLy5 
potentials 27], developed to fit the properties of systems 
with large isospin asymmetry. With account of the ev- 
ident form of FL amplitudes and Eqs. (jnil-ljnil7 one can 
obtain 



?oo = 

?03 = 



p 



2moo 
2mo3 



AiQO, 



M03, 



(14) 
(15) 



6o = (a + 4)^ + 4(q')30' (16) 



£03(P)/03(P) + £33(P)/33(P)}, 



?33 = (c + 4)^ + |(q^)33, 



(17) 
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where the effective nucleon mass moo and effective isovec- 
tor mass mos are defined by the formulae: 



2m, 



00 



A- + ^[3.r+t.(5 + 4..)], (18) 
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ag 
16 



[t2il + 2x2)-ti{l + 2xi)], 



and the renormalized chemical potentials /ioo and /io3 
should be determined from Eqs. 10) • In Eqs. flfil 
and lfT7|l . (q^)3o and (q^}33 are the second order moments 
of the corresponding distribution functions 

(q')30 = ^Eq'/3o(q), (19) 



(q')33 



V 
q 

:^Eq'/33(q). 



(20) 



Thus, with account of the expressions (0) for the dis- 
tribution functions /, we obtain the self-consistent equa- 
tions dl, and |5DJ for the effective chemical 
potentials /i00jM03: FM and AFM spin order parameters 
A£i||, Aqi^, and second order moments (q^)3o, (q^)33- It 
is easy to see, that the self-consistent equations remain 
invariable under a global flip of spins, when neutrons 
(protons) with spin up and spin down are interchanged, 
and under a global flip of isospins, when neutrons and 
protons with the same spin projection are interchanged. 

Let us consider, what differences will be in the case 
of neutron matter. Neutron matter is an infinite nuclear 
system, consisting of nucleons of one species, i.e., neu- 
trons, and, hence, the formalism of one-component Fermi 
liquid should be applied for the description of its proper- 
ties. Formally neutron matter can be considered as the 
limiting case of asymmetric nuclear matter, correspond- 
ing to the isospin asymmetry a = 1. The individual state 
of a neutron is characterized by momentum p and spin 
projection a. The self-consistent equation has the form 
of Eq. ||2J), where all quantities are matrices in the space 
of K = (p, a) variables. The normal distribution function 
and single particle energy can be expanded in the Pauli 
matrices in spin space 



/(P) = /o(p)o-0 + /3(p)cr3, 

£(P) = eo(p)CTo + £3(p)ct3- 



(21) 



The energy functional of neutron matter is character- 
ized by two normal FL amplitudes C/(5'(k) and [/"(k). 
Applying the same procedure, as in Ref. 19], the normal 
FL amplitudes can be found in terms of the Skyrme force 
parameters ti,Xi, (3: 



U^{k)^2to{l-xo) + jg^{l-x^) 

+ ^[<i(l-a;i) + 3*2(1 + a;2)]k2, 

-f -|[<2(l + a;2)-ti(l-xi)]k2 = 
n 



(22) 

(23) 
6„k2. 



With account of Eqs. (|22|l and (|23|) . the normalization 
conditions for the distribution functions can be written 
in the form 
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TT- 



(24) 
(25) 



Here g-^ and g^ are the neutron number densities with 
spin up and spin down and 

h = ]^{n{uj+) + n{uo^)}, c^±=eo±e3, (26) 



Co 
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2m,i 

6 = (an + 5„ — )^^ + — (q )3. 



(27) 
(28) 
(29) 



The effective neutron mass m„ is defined by the formula 
+ |[<i(l-a;i)-f- 3*2(1 + 2:2)], (30) 



g. 



2mn 2mo 

and the quantity (q^)3 in Eq. (|29() is the second order 
moment of the distribution function f^: 



(31) 



Thus, with account of the expressions H26|) and (|27|1 for 
the distribution fimctions /, we obtain the self-consistent 
equations (|24|1 . H25|l . and (|31l) for the effective chemical 
potential /i„, spin order parameter Ag^-^, and second or- 
der moment (q^)3. 



III. PHASE TRANSITIONS IN STRONGLY 
ASYMMETRIC NUCLEAR MATTER 

Early researches on spin polarizability of nuclear mat- 
ter with the Skyrme effective interaction were based on 
the calculation of magnetic susceptibility and finding its 
pole structure 0|, |^ , determining the onset of instabil- 
ity with respect to spin fluctuations. Here we shall find 
directly solutions of the self-consistent equations for the 
FM and AFM spin order parameters as functions of den- 
sity at zero temperature. A special emphasis will be 
laid on the study of strongly asymmetric nuclear matter 
(a ^ 1) while in symmetric nuclear matter FM spin or- 
dering is thermodynamically more preferable than AFM 



If all neutron and proton spins are aligned in one di- 
rection, then for nontrivial solutions of the self-consistent 
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equations we have 



Ag-ff = g, Ag-fi = ag 
(q^)30 ^ 



(32) 



33 



^^gkj.[il + af/' + {l-af/% 
Igkliil + a)'/' - {1 - a)'/% 



where kp = {Sn^g)^^^ is the Fermi momentum of totally 
polarized symmetric nuclear matter. Therefore, for the 
partial number densities of nucleons with spin up and 
spin down one can get 



1 + a 



-Q, Spt 



1-a 



Q, Qpi = Qni = 0. (33) 



If all neutron spins are aligned in one direction and all 
proton spins in the opposite one, then 



A£i|| = ag, Ag^fi = g, 



(34) 



(q')30 
(q')33 



±gkU{l + a)'/' - {1 
^gkUil+a^/' + il 



and, hence. 



-e, Qpi 



-Q, Sp-f = 0ni = 0. (35) 



Now we present the results of numerical solution of 
the self-consistent equations with the effective SLy4 
and SLy5 forces for strongly asymmetric nuclear (a = 
0.95, 0.9, 0.8) and neutron (a = 1) matter. The neutron 
and proton spin polarization parameters n„ and Hp are 
shown in Fig. 1 as functions of density at zero temper- 
ature. Since in totally polarized state the signs of spin 
polarizations are opposite (n„ = l,np = —1), consider- 
ing solutions correspond to the case, when spins of neu- 
trons and protons are aligned in the opposite direction. 
Note that for SLy4 and SLy5 forces, being relevant for 
the description of strongly asymmetric nuclear matter, 
there are no solutions, corresponding to the same direc- 
tion of neutron and proton spins. The reason is that the 
sign of the multiplier t^^—l + 2x3) in the density depen- 
dent term of the FL amplitude Ui, determining spin- 
spin correlations, is positive, and, hence, corresponding 
term increases with the increase of nuclear matter den- 
sity, preventing instability with respect to spin fluctua- 
tions. Contrarily, the density dependent term —t^^g^/S 
in the FL amplitude U3, describing spin-isospin correla- 
tions, is negative, leading to spin instability with the op- 
positely directed spins of neutrons and protons at higher 
densities. 

Another nontrivial feature relates to the density be- 
havior of the spin polarization parameters at large isospin 
asymmetry. As seen from Fig. 1, even small admixture of 
protons leads to the appearance of long tails in the den- 
sity profiles of the neutron spin polarization parameter 
near the transition point to a spin ordered state. As a 



consequence, the spin polarized state is formed much ear- 
lier in density than in pure neutron matter. For example, 
the critical density in neutron matter is g 0.59 fm^'' for 
SLy4 potential and g ~ 0.58 fm~^ for SLy5 potential; in 
asymmetric nuclear matter with a = 0.95 the spin polar- 
ized state arises at p « 0.38 fm~^ for SLy4 and SLy5 po- 
tentials. Hence, even small quantity of protons strongly 
favors spin instability of highly asymmetric nuclear mat- 
ter, leading to the appearance of states with the oppo- 
sitely directed spins of neutrons and protons. As follows 
from Fig. 1, protons become totally spin polarized within 
very narrow density domain (e.g., if a = 0.95, full polar- 
ization occurs at g ~ 0.41 fm~^ for SLy4 force and at 
g « 0.40 fm~^ for SLy5 force) while the threshold densi- 
ties for the appearance and saturation of the neutron spin 
order parameter are substantially different (if a = 0.95, 
neutrons become totally polarized at p « 1.05 fm~^ for 
SLy4 force and at f« 1.02 fm~'^ for SLy5 force). 




1.3 0.4 0.5 0.6 0.7 0.1 



0.9 1.0 1.1 



FIG. 1: Neutron and proton spin polarization parameters as 
functions of density at zero temperature for (a) SLy4 force 
and (b) SLy5 force. 
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FIG. 2: Same as in Fig. but for the second order moments 
{q^)p and (q^)ri, normalized to their values in the totally po- 
larized state. 

Note that the second order moments 

(q2}„ = (q2}„^ - (q2)„^ (36) 

q 
q 

also characterize spin polarization of the neutron and 
proton subsystems. If the solutions (q^)3o and (q^}33 of 
the self-consistent equations are known, then 

(q')n = ^((q')30 + (q')33), 

(q')p-^((q')30-(q')33). 

The values of (q^}n and (q^)p for the totally polarized 
state are 



FIG. 3; Total energy per nucleon, measured from its value in 
the normal state, for the state with the oppositely directed 
spins of neutrons and protons as a function of density at zero 
temperature for (a) SLy4 force and (b) SLy5 force. 



In Fig. El we plot the density dependence of the sec- 
ond order moments (q^)n and (q^)p, normalized to their 
values in the totally polarized state, for different asym- 
metries at zero temperature. These quantities behave 
similar to the spin polarization parameters in Fig.^ i.e., 
there exist long tails in the density profiles of the neutron 
spin order parameter and the proton spin order parame- 
ter is saturated within very narrow density interval. 

To check thermodynamic stability of the spin ordered 
state with the oppositely directed spins of neutrons and 
protons, it is necessary to compare the free energies of 
this state and the normal state. In Fig. 13 the difference 
of the total energies per nucleon of the spin ordered and 
normal states is shown as a function of density at zero 
temperature. One can see that nuclear matter under- 
goes a phase transition to the state with the oppositely 
directed spins of neutrons and protons at some critical 
density, depending on the isospin asymmetry. 
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IV. DISCUSSION AND CONCLUSIONS 

Spin instability is a common feature, associated with 
a large class of Skyrme models, but is not realized in 
more microscopic calculations. The Skyrme interaction 
has been successful in describing nuclei and their excited 
states. In addition, various authors have exploited its 
applicability to describe bulk matter at densities of rel- 
evance to neutron stars The force parameters are 
determined empirically by calculating the ground state 
in the Hartree-Fock approximation and by fitting the 
observed ground state properties of nuclei and nuclear 
matter at the saturation density. In particular, the inter- 
action parameters, describing spin-spin and spin-isospin 
correlations, are constrained from the data on isoscalar 
[2^133 and isovector (giant Gamow-Teller) |31j j52] spin- 
flip resonances. 

In a microscopic approach, one starts with the bare 
interaction and obtains an effective particle-hole inter- 
action by solving iteratively the Bethe-Goldstone equa- 
tion. In contrast to the Skyrme models, calculations with 
realistic NN potentials predict more repulsive total en- 
ergy per particle for a polarized state comparing to the 
nonpolarized one for all relevant densities, and, hence, 
give no indication of a phase transition to spin ordered 
state. It must be emphasized that the interaction in the 
spin- and isospin-dependent channels is a crucial ingredi- 
ent in calculating spin properties of isospin asymmetric 
nuclear matter and different behavior at high densities 
of the interaction amplitudes, describing spin-spin and 
spin-isospin correlations, lays behind this divergence in 
calculations with the effective and realistic potentials. 

In this study as a potential of NN interaction we choose 
SLy4 and SLy5 Skyrme effective forces, which were con- 
strained originally to reproduce the results of micro- 
scopic neutron matter calculations (pressure versus den- 
sity curve) [231 . Besides, in the recent publication l28| 
it was shown that the density dependence of the nuclear 
symmetry energy, calculated up to densities g 'SSqq with 
SLy4 and SLy5 effective forces (together with some other 
sets of parameters among the total 87 Skyrme force pa- 
rameterizations checked) gives the neutron star models 
in a broad agreement with the observables, such as the 
minimum rotation period, gravitational mass-radius re- 
lation, the binding energy, released in supernova collapse, 
etc. This is important check for using these parameter- 
izations in high density region of strongly asymmetric 
nuclear matter. However, it is necessary to note, that 
the spin-dependent part of the Skyrme interaction at 
densities of relevance to neutron stars still remains to 



be constrained. Probably, these constraints will be ob- 
tained from the data on the time decay of magnetic field 
of isolated neutron stars :33| . In spite of this shortcom- 
ing, SLy4 and SLy5 effective forces hold one of the most 
competing Skyrme parameterizations at present time for 
description of isospin asymmetric nuclear matter at high 
density while a Fermi liquid approach with Skyrme effec- 
tive forces provides a consistent and transparent frame- 
work for studying spin instabilities in a nucleon system. 

In summary, we have considered the possibility of 
phase transitions into spin ordered states of strongly 
asymmetric nuclear matter within the Fermi liquid for- 
malism, where NN interaction is described by the Skyrme 
effective forces (SLy4 and SLy5 potentials). In contrast 
to the previous considerations, where the possibility of 
formation of FM spin polarized states was studied on 
the base of calculation of magnetic susceptibility, we ob- 
tain the self-consistent equations for the FM and AFM 
spin polarization parameters and solve them in the case 
of zero temperature. It has been shown in the model with 
SLy4 and SLy5 effective forces, that strongly asymmet- 
ric nuclear matter undergoes a phase transition to the 
spin polarized state with the oppositely directed spins 
of neutrons and protons, while the state with the same 
direction of the neutron and proton spins does not ap- 
pear. An important peculiarity of this phase transition 
is the existence of long tails in the density profile of the 
neutron spin polarization parameter near the transition 
point. This means, that even small admixture of protons 
to neutron matter leads to the considerable shift of the 
critical density of spin instability in the direction of low 
densities. In the model with SLy4 effective interaction 
this displacement is from the critical density g ~ 3.7£'o 
for neutron matter to g w 2.4£»o for asymmetric nuclear 
matter at the isospin asymmetry a = 0.95, i.e. for 2.5% 
of protons only. As a result, the state with the oppo- 
sitely directed spins of neutrons and protons appears, 
where protons become totally polarized in a very narrow 
density domain. This picture is different from the case 
of symmetric nuclear matter, where the FM spin config- 
uration is thermodynamically more preferable, than the 
AFM one j^J. Obtained results may be of importance 
for the description of thermal and magnetic evolution of 
pulsars, whose core represents strongly asymmetric nu- 
clear matter. 
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